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Abstract 

We study initial-boundary value problems for linear evolution equations of arbitrary spa- 
tial order, subject to arbitrary linear boundary conditions and posed on a rectangular 1-space, 
1-time domain. We give a new characterisation of the boundary conditions that specify well- 
posed problems using Fokas' transform method. We also give a sufficient condition guarantee- 
ing that the solution can be represented using a series. 

The relevant condition, the analyticity at infinity of certain meromorphic functions within 
particular sectors, is significantly more concrete and easier to test than the previous criterion, 
based on the existence of admissible functions. 

1 Introduction 

In this work, we consider 

The initial-boundary value problem Il{n, A, a,h, qo): Find q £ C°°([0, 1] x [0, T]) which 
satisfies the linear, evolution, constant-coefficient partial differential equation 

dtqix,t) + a{-id,)"qix,t) =0 (1.1) 

subject to the initial condition 

q{x,0)=qo(x) (1.2) 

and the boundary conditions 

A {dr\{o, t), ar t), a^-'g(o, t), ar 'g(i, t), • • • , g(o, t), g(i, t)Y = h{t), (1.3) 

where the pentuple (n. A, a, h, qo) £ N x R"x2n ^ ^ (C°°[0, T])" x C°°[0, 1] is such that 
(ni) the order n ^ 2, 

(112) the boundary coefficient matrix A is in reduced row-echelon form, 

(113) if n is odd then the direction coefficient a = ±i, if n is even then a = e'* for some 

ee[-7v/2,7T/2], 

(114) the boundary data h and the initial datum qo are compatible in the sense that 

A(g("-^)(0),g("-^'(l),,,r''(0),?^"-''(l),...,9o(0),go(l))^=MO). (1-4) 

Provided 11 is well-posed, in the sense of admitting a unique, smooth solution, its solution 
may be found using Fokas' unified transform method [5l[7]. The representation thus obtained 
is a contour integral of transforms of the initial and boundary data. Certain problems, for 
example those with periodic boundary conditions, may be solved using classical methods such 
as Fourier's separation of variables JQ], to yield a representation of the solution as a discrete 
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Fourier series. By the well-posedness of H, these are two different representations of the same 
solution. 

For individual examples, Pelloni 13 and Chilton [5] discuss a method of recovering a series 
representation from the integral representation through a contour deformation and a residue 
calculation. Particular examples have been identified of well-posed problems for which this 
deformation fails but there is no systematic method of determining its applicability. 

Pelloni [12] uses Fokas' method to decide well-posedness of a class of problems with un- 
coupled, non-Robin boundary conditions giving an explicit condition, the number that must 
be specified at each end of the space interval, whose validity may be ascertained immediately. 
However there exist no criteria for well-posedness that are at once more general than Pelloni's 
and simpler to check than the technical 'admissible set' characterisation of [7]. 

The principal result of this work is a new characterisation of well-posedness. The condition 
is the decay of particular integrands within certain sectors of the complex plane. Indeed, let 
D = {p G C : Re(ap") < 0}. Then 

Theorem 1.1. The problem Il{n, A,a, h, qo) is well-posed if and only if rij{p) is entire and 
the ratio 



vAp) 



as p CO from within D, away from the zeros of Apde . 



(1.5) 



for each j. 

We provide a small contribution to Fokas' method, making it fully algorithmic. We express 
the solution in terms of the PDE characteristic determinant, Apde , the determinant of the 
matrix 



Aj(p) = <^ 



reJ + 



-He" 



)/2" 



^(„_l_r)(fc_l)^^^^(J^-_l)/2-r 



rS J- 



(1.6) 



V,/2(P) ^-w'' 
re./ + 



/2" 



_(„_l_r)(fc-l)^^^-^Jj./2-r 



Jj even. 



rS J- 



The matrix A appears in the generalised spectral Dirichlet to Neumann map derived in Sec- 
tion [21 The application of the map to the formal result Theorem 12.11 yields the following 
implicit equation for q, the solution of 11. 

Theorem 1.2. Let Tl{n,A,a,h,qo) be well-posed with solution q. Then q(x,t) may be ex- 
pressed in terms of contour integrals of transforms of the boundary data, initial datum and 
solution at final time as follows: 



2nq{x, t) 



I'' 



px — ap t 



0(p)-e"'' ^% (p) 



dp 



an- 



E 



ApDE (p) 

I — ( 

ApDE (p) 



dp, (1.7) 



where the sectors = DnC^ and D = {p G C : Re(ap") < 0}. 
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Equation (II. 7p gives only an implicit representation of the solution as the functions rjj are 
defined in terms of the Fourier transform of the solution evaluated at final time, which is not 
a datum of the problem. Nevertheless the importance of the PDE characteristic determinant 
is clear. The integrands are meromorphic functions so q depends upon their behaviour as 
p ^ 00 from within and upon their poles, which can only arise at zeros of Apde • It is 
the behaviour at infinity that is used to characterise well-posedness in Theorem 11.11 the proof 
of which is given in Section O 

In Section |4] we derive two representations of the solution of an initial-boundary value 
problem. Let {(Jk)keTi be a sequence containing each nonzero zero of Apde precisely once and 
define the index sets 

= {fc G N : CTfe G R}, 
{fc G N : ImcTfc > 0}, 
{fc G N : Imafc < 0}. 

Then the following theorems give representations of the solution to the problem 11. 

Theorem 1.3. Let the problem Il{n,A,a,h,qo) be well-posed. Then the solution q may be 
expressed using contour integrals of transforms of the initial and boundary data by 

2 "^^^ ,tj^ 

+ 2 ^f, Apde (p) ^ J ^J,-' ^ Apde (p) 

" ^ {.Sir//.} 

Theorem 1.4. Let a = ±i and let the problems 11 — n(n. A, a, h, go) and 11' = n(n, A, —a, h, go) 
be well-posed. Then the solution q of II may be expressed as a discrete series of transforms of 
the initial and boundary data by 

^ ^ipx — ap^t 

2 p=<yk Apde (p) 

+ 9 E -1^ — ur E '^j(^) 

2, ^ — ' P=o-fc ApDElp) ^ — ' 

The final integral term in both equations (|1.8p and (|1.9p depends upon H, a linear com- 
bination of t-transforms of the boundary data which evaluates to if /i = 0. Hence if 11 is 
a homogeneous initial-boundary value problem then the final term makes no contribution to 
equations (|1.8|l and (|1.9|l . 

Special cases of Theorem 11.31 have appeared before but the representations differ from 
equation (|1.8|l . The result is shown for several specific examples in [51 [13] , including a second 
order problem with Robin boundary conditions. For simple boundary conditions, the result is 
mentioned in Remark 6 of 7 and Lemma 4.2 of [T^] contains the essence of the proof. Unlike 
earlier forms, equation (|1.8p represents q using discrete series as far as possible; only the parts 
of the integral terms that cannot be represented as series remain. This may not have any 
advantage for computation but is done to highlight the contrast with equation (|1.9p . 

In Theorem 11.41 the well-posedness of EL' is used to show that the first two integral terms 
of equation (|1.8p evaluate to zero. Under the map a i— >■ —a, D maps to E, the interior of its 
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complement; we exploit this fact together with Theorem 1 1.1 1 to show the decay of 



Up) 



as p — )■ oo from within E. 



ApDE (p) 



This maximally generalises of the arguments of Pelloni and Chilton in the sense that the 
deformation of contours cannot yield a series representation of the solution to 11 if H' is ill- 
posed. 

Theorem 1 1.1 1 is useful because it reduces the complexity of the analysis necessary to prove 
that a particular initial-boundary value problem is well-posed but its use still requires some 
asymptotic analysis. It would be preferable to give a condition that may be validated by 
inspection of the boundary coefficient matrix and is sufficient for well-posedness. We discuss 
such criteria in Section O 

Section [S] also contains a proof of the following result, complementing Theorem 11.41 This 
theorem highlights the essential difference between odd order problems, whose well-posedness 
depends upon the direction coefficient, and even order problems, whose well-posedness is 
determined by the boundary coefficient matrix only. 

Theorem 1.5. Let n be even and a — ±i. Using the notation of Theorem \1.4\ the problem 
n is well-posed if and only ifH is well-posed. 

In Section [6] we investigate the PDE discrete spectrum, the set of zeros of the PDE charac- 
teristic determinant. We prove a technical lemma, describing the distribution of the at which 
is used in the earlier sections. Under certain conditions we are able to exploit symmetry argu- 
ments to improve upon the general results Langer presents [IT] for the particular exponential 
polynomials of interest. 



In Section [27T] we give the standard results of Fokas' unified transform method in the notation 
of this work. In Section [2.21 we state and prove Lemma [2.61 the generalised spectral Dirichlct 
to Neumann map. In Section 12.31 we apply the map to the formal results of Section 12.11 
concluding the proof of Theorem 11.21 The latter two sections contain formal definitions of 
many of the terms and much of the notation used throughout this work. 

2.1 Fokas' method 

The first steps of Fokas' transform method yield a formal representation for the solution of 
the initial-boundary value problem, given in the following 

Theorem 2.1. Let the initial-boundary value problem Il(n, A, a, h, qo) be well-posed. Then its 
solution q may be expressed formally as the sum of three contour integrals, 



2 Implicit solution of IBVP 





where 




e 



ap s 




hit) 



diq{0,t). 



9 At) 



diqiht), 



(2.2) 





-0+1) 
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The above theorem is well established and its proof, via Lax pair and Riemann-Hilbert 
formalism, appears in (4] [5l [?]■ We state it here without proof to highlight the difference in 
notation to previous publications. We use p to denote the spectral parameter, in place of k 
in the earlier work. We use /, and Qj exclusively to denote the boundary functions; even for 
simple boundary conditions in which some of the boundary functions are equal to boundary 
data we denote the boundary data separately by ■ 

The transformed boundary functions are the 2n unknowns in equation (|2.ip . of which 
at most n may be explicitly specified by the boundary conditions H1.3|) . To determine the 
remaining n or more we require a generalised Dirichlet to Neumann map in the form of 
Lemma 12.61 This is derived from the boundary conditions and the global relation. 

Lemma 2.2 (Global relation). Let Il(n,A,a,h,qo) be well-posed with solution q. Let 

qT{p)= f e-''\{y,T)Ay 
Jo 

be the usual spatial Fourier transform of the solution evaluated at final time. Then the trans- 
formed functions qo, qr, fj and gj satisfy 

n-l 

E ^^(P^ [f^^p) - ^"^'Mp)) = Up) - e^'"^qT{p), p e c. (2.3) 

The global relation is derived using an application of Green's Theorem to the domain 
[0, 1] X [0, T] in the aforementioned publications. As the f-transform, 

X{p)= r e'^''^'X{t)dt, (2.4) 
Jo 

is invariant under the map p i— >■ exp {2jTii/n)p for any integer j, the global relation provides a 
system of n equations in the transformed functions to complement the boundary conditions. 

2.2 Generalised spectral Dirichlet to Neumann map 

We give a classification of boundary conditions and formally state the generalised spectral 
Dirichlet to Neumann map. 

Notation 2.3. Consider the problem \l{n, A,a,h, qo), which need not be well-posed. Define 
uj = exp {2ni/n). Define the boundary coeffictents Okj, Pkj to be the entries of A such that 

/ ai „_i 

a2n-l 

Tin— 1 Pnn—1 — 2 Pn n — 2 ■■■ an Pn 0/ 

We define the following index sets and functions. 

J+ = {j £ {0, 1, . . . , n — 1} such that akj is a pivot in A for some k}, the set of columns 
of A relating to the left of the space interval which contain a pivot. 

J~ = {j £ {0, 1, ... ,n — 1} such that Pkj is a pivot in A for some k}, the set of columns 
of A relating to the right of the space interval which contain a pivot. 

= {0, l,...,n— 1}\J"'", the set of columns of A relating to the left of the space interval 
which do not contain a pivot. 

J~ = {0, 1, . . . ,n — 1} \ J~ , the set of columns of A relating to the right of the space 
interval which do not contain a pivot. 

J = {2j + 1 such that j £ J^} U {2j such that j £ J }, an index set for the boundary 
functions whose corresponding columns in A do not contain a pivot. Also, the decreasing 
sequence {JjYj^i of elements of J. 

J' = {2j + 1 such that j £ J+} U {2j such that j £ J~} ^ {0,1, . . . ,2n - 1} \ J, an index 
set for the boundary functions whose corresponding columns in A contain a pivot. Also, the 
decreasing sequence (Jj)"^! of elements of J'. 



Pm- 

P2n- 



ai , 

a2r 



Pi n-2 
P2 n-2 



aio Pio\ 
0.2 P2 



A. 



(2.5) 
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The functions 

V{p) = (Vi(p), V2(p), . . . , VAp) = 



/(j,-i)/2(p) >/j odd, 
gj,/2{p) Jj even, 



the boundary functions whose corresponding columns in A do not contain a pivot. 
The functions 



W(p) = (Wiip), W2{p), W„{p)f, WAp) = 



/(j'-i)/2(p) -/j' odd, 
9j'/2{p) Jj even. 



the boundary functions whose corresponding columns in A contain a pivot. 
{Jj')j^j+: a- sequence such that ciy+ ^ is a pivot in A when j £ J+. 

{ Jj")j^ J- , a sequence such that Pj- ^ is a pivot in A when j £ J~ . 

Definition 2.4 (Classification of boundary conditions). The boundary conditions of the prob- 
lem n(n. A, a, h, qo) are said to be 

1. homogeneous if h — Q. Other-wise the boundary conditions are inhomogeneous. 

2. uncoupled if 

if akj is a pivot in A then r = V r and 
if Pk j is a pivot in A then = V r. 

Otherwise we say that the boundary conditions are coupled. 

3. non-Robin if 

V fc G {1,2, . . . ,n}, ifokjT^OorPkj^O then Uk r = Pkr ^ W r ^ j, 

that IS each contains only one order of partial derivative. Otherwise we say that boundary 
condition is of Robin type. Note that whether boundary conditions are of Robtn type or 
not is independent of whether they are coupled, unlike Duff's definition 

4- simple if they are uncoupled and non-Robm. 

The terms 'generalised' and 'spectral' are prefixed to the name 'Dirichlet to Neumann map' 
of the Lemma below to avoid confusion regarding its function. 

Generalised: The boundary conditions we study are considerably more complex than 
those considered in piHirn[51 ll2lll3| . Indeed, as A may specify any linear boundary conditions, 
the known boundary functions may not be 'Dirichlet' (zero order) and the unknown boundary 
functions need not be 'Neumann' (first order). Further, if A has more than n non-zero entries 
then the lemma must be capable of expressing more than n unknown boundary functions in 
terms of fewer than n known boundary data. 

Spectral: Owing to the form of equation (|2.1|l we are interested not in the boundary 
functions themselves but in their t-transforms, as defined in equations (|2.2|l . It is possible, 
though unnecessarily complicated, to perform a generalized Dirichlet to Neumann map in real 
time and subsequently transform to spectral time but, as the global relation is in spectral time, 
to do so requires the use of an inverse spectral transform. Instead, we exploit the linearity of 
the t-transform (|2.4p . applying it to the boundary conditions, and derive the map in spectral 
time. 

The crucial component of the lemma is given in the following 

Definition 2.5. Let n(n. A, a, h, qo) be an inittal-boundary value problem having the properties 
(ni)-(n4) but not necessarily well-posed. We define the PDE characteristic matrix A{p) by 
equation p.6p and the PDE characteristic determinant to be the entire function 

Apde(p) = det^(p). (2.6) 

Lemma 2.6 (Generalised spectral Dirichlet to Neumann map). Let n(n, A, a, h,, go) be well- 
posed with solution q. Then 
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1. The vector V of transformed boundary functions satisfies the reduced global relation 
A{p)V{p) = U{p) - e"""^ 



(2.7) 



where 

U{p) = {u{p,l),u{p,2),...,u{p,n)f, (2.8) 
u{p,k) = qo{uj^'^p) - ^ c;(cj''"V)^j+(p) + e"" " ^ ci{u)''~^ p)hj-{p) (2.9) 

and /ij is t/ie function obtained by applying the t-transform p.4|l to the boundary datum 
hj. 

2. The PDE characteristic matrix is of full rank, is independent of h and go o,nd differing 
values of a only scale A by a nonzero constant factor. 

3. The vectors V and W of transformed boundary functions satisfy the reduced boundary 
conditions 

W{p)^[Mp),h2{p),---X{p)y -AV{p), (2.10) 
where the reduced boundary coefficient matrix is given by 

f«.U.-i)/2 J^odd, ^2.11) 
Jj even. 

Proof. Applying the t-transform (12. 4p to each line of the boundary conditions (|1.3p yields a 
system of n equations in the transformed boundary functions. As A is in reduced row-echelon 
form it is possible to split the vector containing all of the transformed boundary functions into 
the two vectors V and W , justifying the reduced boundary conditions. 
The reduced boundary conditions may also be written 

fj (p) = h j+ (p) - ^ a j+ Jr (p) ~ J2 h+ r9r (p) , for j G J+ and (2.12) 

re. 7+ rGJ- 

9J ip) = hj^ (p) - J2 «J- Mp) - E I^J- Mp), for j e J- . (2.13) 

reJ+ rej- 

As the t-transform is invariant under the map p i-^ ui-' p for any integer j, the global relation 
Lemma 12.21 vields the system 

n—l n— 1 

X] Cj(p)t^'""^"-''''/5(p) - e""^''''cj(p)a;<""^"^''pj(p) = go('^» - e'"'"^gT(w''p), 
j=o i=o 

for r G {0, 1, . . . , n - 1}. Using the fact J+ U J+ = J " U J ~ = {0, 1, . . . , n - 1} we split the 
sums on the left hand side to give 

c,{p)j-'-^^^f,{p)+ J2 cApW--'-'^^Jap) 

je.J+ 

- J2 e-'-''''c,(p)c.("-^-^>5.(p)- J2 e-'-^'cApW-'-'^'^Mp) 
js./- jeJ- 

= go(t^ p)-e qriuJ p), 

for r G {0,1,..., n — 1}. Substituting equations (I2.12|l and (|2.13|l and interchanging the 
summations we obtain the reduced global relation. 

The latter statement of (ii) is a trivial observation from the form of the PDE characteristic 
matrix. A full proof that A is full rank is given in the proof of Lemma 2.17 of [II]. □ 
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2.3 Applying the map 

We solve the system of linear equations H2.7[l for V using Cramer's rule hence, by equa- 
tion (|2.10|) . determining W also. 

Notation 2.7. Denote by Q{p) the determinant of the matrix obtained by replacing the 
j"^ column of the PDE characteristic matrix with the vector U{p) and denote by fjj{p) the 
determinant of the matrix obtained by replacing the j"" column of the PDE characteristic 
matrix with the vector (gT(p), <?t(i^p), ■ • ■ , qT{(^"~^ p))'^ for j e {1, 2, . . . , n} and p £ C. Define 

n 

Zip) ^hj-n{p) - ^Aj-„kCkip), 



%{p) = hj-^{p) - ^ Aj_„fc?7fe(p), 



(2.14) 



for j £ {n + l,n + 2, . . . , 2n} and p £ C. Define 



0(p) 



C(J,-l)/2(p)Cj(p) 

CJ,/2(P)C;(P) 

C(J'_^-l)/2(p)Cj(p) 

Cj' ^/2(P)C;(P) 



C(j,~i)/2[p)riAp) 

Cj^/2{p)rjAp) J] 

%J'^_^-i)/2{p)m{p) J'o 

CJ'_„/2(P)%(P) 



Jj odd, 



odd. 



(2.15) 



J'j-n even. 



for p £ C and define the index sets 

J"^ = {j ■ Jj odd} U{n + j: jj odd}, 
J~ = {j : Jj even} U {n + j : Jj even}. 

The generalised spectral Dirichlet to Neumann map Lemma 
expressions for the transformed boundary functions: 



and Cramer's rule yield 



Ci(p) 



vAp) 



:(P) 



C(jj-1)/2(P)/(J,-1)/2(P) 
Cj^/2{p)gjj/2{p) Jj 
C(Jj_^-l)/2(p)/(Jj_„-l)/2(p) Jj 
Cj'_ /2(p)gj'_ /2(P) 



Jj odd. 



7' 

''j—n 



odd, 
even, 



(2.16) 



hence 



j=o jeJ+ 

n-l 

^<^3{p)gj{p) = E 

.7 = jSJ- 



G(P)- 



%(P) 



;(P) 



0(p) 



%(P) 



:(P) 



Substituting these equations into Theorem 12.11 completes the proof of Theorem 11.21 

Remark 2.8. There are several simplifications of the above definitions for specific types of 

boundary conditions. 

If the boundary conditions are simple, as studied in 12,, then A = Q. Hence, if the 
boundary conditions are simple and homogeneous then C,j — rjj = for each j > n. 

Non-Robin boundary conditions admit a significantly simplified form of the PDE charac- 
teristic matrix; see equation (2.2.5) of [14) . 

For homogeneous boundary conditions, rjj is Q with qr replacing go- 
Remark 2.9. It is possible to extend the results above to initial-boundary value problems for 
a more general linear, constant-coefficient evolution equation. 



dtq{x,t) + '^aj{-idxy q{x,t) = 0, 

j=0 



(2.17) 
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with leading coefficient a,i liaving the properties of a. In this case the spectral transforms 
must be redefined with X]"^o "'iP'' replacing ap" and the form of the boundary coefficient 
matrix also changes. The w'^ appearing in equation p.6p represent a rotation by 2X7r/n, 
corresponding to a map between simply connected components of D. The partial differential 
equation (|2.17|l has dispersion relation X^"^o "jP" ^ ^'^^ simply a union of sectors but a 
union of sets that are asymptotically sectors; see Lemma 1.1 of [9]. Hence we replace cj'^ with 
a biholomorphic map between the components of D. 



3 New characterisation of well-posedness 

This section provides a proof of Theorem 11.11 The first subsection justifies that the decay 
condition is satisfied by all well-posed problems. The second subsection proves that the decay 
condition is sufficient for well-posedness. 

We clarify the definitions of D and E from Section [T] By Lemma 16.11 there exists some 
e > such that the pairwise intersection of closed discs of radius e centred at zeros of Apde 
is empty. We define 

D = D\yjB{cjk,e), E = E\y]B{ok,e). 



3.1 Well-posedness decay 

As the problem is well-posed, the solution evaluated at final time qt £ C°°[0, 1] hence qx and 
Tjj are entire. Similarly, fk,gk € C°°[0,r] hence fk,gk are entire and decay as p — >■ oo from 
within D. Hence, by equation (|2.16p . 

^-y-f;>i^) (3-1) 

ApDE(p)Cfc(p) 

is entire and decays as p — >■ from within D for each j G {1, 2, . . . , 2n}, where k depends upon 
j- 

We define the new complex set 

V = {peD such that - Re(ap"T) > 2n|p!}. 



As C -D, the ratio (|3.ip is analytic on O and decays as p — > oo from within T>. For 
p £ {1, 2, . . . , n}, let Dp be the p"^ simply connected component of D encountered when moving 
anticlockwise from the positive real axis and let Dp = Dn Dp. Then for each p £ {l,2,...,n} 
there exists R > such that the set 

^DpHV^ \B{0,R) 

is simply connected, open and unbounded. 

By definition, Apde (p) is an exponential polynomial whose terms are each 

W(p)e-'S«e^"'''' 

where is a monomial of degree at least 1 and Y C {0,l,2,...,n— 1} is an index set. Hence 

o(e"''''p~^) as p — > oo or as p — )■ 0. 



ApDE (p) 

As Q and r/j also grow no faster than o(e"''''), the ratios 

— o(e ' 'p ), as p — >■ oo. 



CAp) Vjip) _ 2n|p| 



ApDE (p)Cfc (p) ' ApDE (p)cfc (p) 

Hence the ratio 



7^ (3.2) 
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decays as p — >■ oo from within V and away from the zeros of Apde ■ However the ratio 

0(p) 

ApDE {p)Ck{p) 



(3.3) 



is the sum of ratios p.l|) and 1)3. 2[1 hence it also decays as p — > oo from within D' and away 
from the zeros of Apde . 

The terms in each of C,j (p) and Apde (p) are exponentials, each of which either decays or 
grows as p — > oo from within one of the simply connected components Dp of D. Hence as 
p — > oo from within a particular component Dp the ratio 1)3. 3|) either decays or grows. But, as 
observed above, these ratios all decay as p — >■ oo from within each T>p. Hence the ratio 1)3. 3|) 
decays as p — >■ oo from within Dp. 

Now it is a simple observation that the ratio 

(3.4) 



ApDE {p)Ck{p) 

must also decay as p — >■ oo. Indeed ratio (|3.4p is the same as ratio (|3.3p but with (Jt('^'^~^p) 
replacing u{p,k) and, as observed above, qt G C°°[0, 1] also. Finally, the exponentials in rjj 
and ApDE ensure that the ratio 

vAp) 



ApDE (p) 



(3.5) 



also decays as p — >■ oo from within Dp. Indeed the transforms that multiply each term in Tjj 
ensure that the decay of ratio (|3.4|) must come from the decay of ratio (|3.5|) . not from l/cfc(p). 

3.2 Decay =^ well-posedness 

Many of the definitions of Section [2] require the problem n(n. A, a, h, qo) to be well-posed. The 
statement of the following Lemma clarifies what is meant by rjj when 11 is not known to be 
well-posed a priori and the result is the principal tool in the proof of Theorem 11.11 

Lemma 3.1. Constder the problem 11 (n, A, a, h, qo) with associated PDF characteristic matrix 
A whose determinant is Apde • Let the polynomials Cj be defined by Cj{p) = — ap"(ip)^'-'^^' . 



Let U : C — > C 6e defined by equation (|2.8|) and let A G R"*^" be defined by equation H2.10|l . 
Let Q,rij : C ^ C be defined by Notation \2. 7[ where qr : [0, 1] — > C is some function such 
that Tjj is entire and the decay condition H1.5|) is satisfied. Let the functions fj ,gj : C — >■ C be 



defined by equation ()2.16|) . Let fj , Qj : [0, T] ^ C be the functions for which 

IAp)- r e'^''"\f,{t)dt, g,{p)^ r e'"'"'g,it)dt, p G C. (3.6) 
Jo Jo 

Then {fj,gj : j G {0, 1, . . . ,n — 1}} is an admissible set in the sense of Definition 1.3 of J^. 



Proof. By equation (|2.16p and the definition of the index sets J in Notation 12.71 we may 
write equations (1.13) and (1.14) of [71 as 

Hp)-J: '^^'^-'^':Y'^ , (3.7) 
ApDE (p) 



jeJ+ 

Gip)^j:^-M^^^^. (3.8) 

ApDE (p) 

By Cramer's rule and the calculations in the proof of Lemma 12.61 equation (1.17) of ;7 is 
satisfied. 

As rjj is entire, qr is entire so, by the standard results on the inverse Fourier transform, 
qr : [0, 1] — >■ C, defined by 

1 



<1t{x) = ^ / e""'qT(p)dp, 
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is a C°° smooth function. 

We know (j is entire by construction and rjj is entire by assumption hence F and G are 
meromorphic on C and analytic on D. By the definition of D and the decay assumption 

' as p -> oo from within D. 



ApDE (p) 

As go and hj are entire so is U. As qr is also entire and the definitions of Q and rjj differ only 
by which of these functions appears, the ratio ("fc(p)/ApDE (p) — >■ as p — > oo from within D 
also. This establishes that 

0(p) -e''P^r?,(p) _^ as p ^ oo from within D. 
ApDE (p) 

Hence, by equations (|3.7p and (|3.8|l . F{p), G{p) — >■ as p — )■ oo within _D. 
An argument similar to that in Example 7.4.6 of [Ij yields 



gAt) = ^f f P^e""''"*G(p)dp. 
J an 



Because F(p),G{p) — >■ as p — !> cxj within D, these definitions guarantee that fj and are 
C°° smooth. 

The compatibility of the fj and with go is ensured by the compatibility condition 
(n4). □ 

The desired result is now a restatement of Theorems 1.1 and 1.2 of [?]• For this reason we 
refer the reader to the proof presented in Section 4 of that publication. The only difi'erence is 
that we make use of Lemma |3. II in place of Proposition 4.1. 



4 Representations of the solution 

The proofs of Theorems 11.31 and 11.41 are similar calculations. In Section 14.11 we present the 
derivation of the series representation and, in Section [4.21 note the way this argument may be 
adapted to yield the integral representation. We derive the result in the case n odd, a — i; 
the other cases are almost identical. 

4.1 Series Representation 

As n is well-posed. Theorem 11.21 holds. We split the latter two integrals of equation (|1.7|l into 
parts whose integrands contain the data, that is Q, and parts whose integrands contain the 
solution evaluated at final time, that is rjj. 

f g,p.+,p"(T-o y vAp) dp + I / + f\ e»''(-i)-p"' y J^^P] ^ dp 
Jbd+ .f7iApDE(p) \JaE- Jm ^ Apde (p) 

+ / g«P(-i)+«P"(T-t) y- J)Ap) 

ho- ApDE (p) 

As n' is well-posed. Theorem 11.11 ensures the ratios 

■ ^-'^^^ ^ — > as p — >■ 00 from within Z)', 

Ai>DE (P) 
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for each j. By definition E — D' and, by statement (ii) of Lemma 12.61 the zeros of A' are 
precisely the zeros of A hence E — D' . Define Cj (p) to be the function obtained by replacing 
q'j'{uj''~^ p) with u{p,k) in the definition of ri'j{p). As q'^, qo and hj are all smooth functions, 
5j has precisely the same decay properties of ri'j. But = by definition. Hence the 
well-posedness of fl' is equivalent to 

^ as p -> oo from within E, (4.2) 



ApDE (p) 

for each j. The decay property obtained by applying Theorem II. fl directly to fl together 
with the decay property (|4.2p permits the use of Jordan's Lemma to deform the contours of 
integration over and i?^ in equation (|4.1|l to obtain 



27rqix,t) = / e'''^-"'"*go(p)dp+ \ [ - [] e'"^"'^"* ^ 



^dp 



ApDE (p) 

je j+ 

Jd{D + \D+) j^+^P'O'^^P^ 

+ W.,._,._, + /.h"^'^-^"'"'E^dp 

ip{x-l)+ip"{T-t) % (P) 



a(B-\E-) 



+ / g.p(.-i)+.p"(T-t) (4 3) 

Jd{D-\D-) Apde(p) 

Indeed, (^j , r^j and Apde are entire functions hence the ratios can have poles only at the zeros 
of ApDE , neighbourhoods of which are excluded from Z)^ and E^ by definition. Finally, the 
exponential functions in the integrands each decay as p — )■ oo from within the sectors enclosed 
by their respective contour of integration. 

The right hand side of equation (|4.3p is the sum of three integrals over R and four others. 
The former may be combined into a single integral using the following lemma, whose proof 
appears at the end of this section. 

Lemma 4.1. Let Il{n, A, a, qo, h) be well-posed. Then 



E Cj{p)~e-^^ E 0(p) = Apde(p) 

jeJ+ jeJ- 



(4.4) 



where 

H{p)= E c,(p)fty+(p)-e-'^ E C.(P)V(P). 

jeJ- 

The other integrals in equation 1)4. 3|) are around the boundaries of discs and circular sectors 
centred at each zero of Apde • Over the next paragraphs we combine and simplify these 
integrals to the desired form. 

Consider a £ such that Apde (o") = 0. Then the fourth integral on the right hand side 
of equation (|4.3|) includes 



C(a,e) Apde (p) }c(.,e) Apde (p) 

the equality being justified by the following lemma, whose proof appears at the end of the 
section. 

Lemma 4.2. LetX\[n, A,a,qo,h) be well-posed. Then for each j £ {I, 2, . . . , 2n}, the functions 

0(p)-e°""^^.('^) V- 0(p)-e''^"^%(a) 

Apde(p) ' ^ Apde(p) ^ • ' 

are entire. 
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Consider a G {dD) n C+ such that Apde (cr) = 0. Define = d{B{a,e) n D) and = 
d{B{a,e) H -E). Then the second and fourth integrals on the right hand side of equation H4.3|l 
include 



/ -T Cj(p)dp and 

ipB ApDE (p) 



/pB ApDE (p) 

ApDE (p) JpdApde(P) 

j e j+ jeJ+ 
respectively, by Lemma [4.21 The sum of the above expressions is 

r ^ipx-ip"t 

J a.,.) ApDE (p) 

Consider / cr £ R such that Apde (cr) = 0. Define r° = d{B{a,e) n D) and let 
= d{B{a, e)nE). Then the fourth and fifth integrals on the right hand side of equation 
include 



/ e-+-"(--)E A^dp=/ ^:!ZQ^C.(p)dp and 
Jro ApDE (p) JrD ApDE (p) , 

JpE ApDE (p) .-^ 
respectively, by analyticity and Lemma 14.21 The sum of the above expressions is 

r ipx — ip^t r ipx — ip^t I \ 

Apde(p) Apde(p) y 

Similar calculations may be performed for cr G iJ^, D~, H C^, {0}. Define the index 
set A'* C N by fc G J^* if and only if crfc G R. For each k £ K^' define Ffe = d{B{ak,e) D C"). 
Then, substituting the calculations above and applying Lemma [4.11 equation (|4.3p yields 

^ ^ r ipx-ip"t ^ 



ip{x — l) —ip^t 

+ E L . A„..r.^ EO(p)dp 



C{<Tfc,£) Apde(p) 



A/ VApde(p) i 

A residue calculation at each at completes the proof. 



4.2 Integral Representation 

As n is well-posed, equation l|4.1|) holds but, as Il{n, A, —a, h, qo) may not be well-posed, it 
is not possible to use Jordan's Lemma to deform the second and fifth integrals on the right 
hand side over E. However it is still possible to deform the fourth and seventh integrals over 
D. Hence two additional terms appear in equation (|4.3p . 

f ^^P--^p-ty ^Ap) , f ^^pix-l)-^p-t y Cjip) . 

Jqe+ f-; Apde (p) J ais- t1 ^pde (p) 

The remainder of the derivation is unchanged from that presented in Section [4.11 
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4.3 Proofs of technical lemmata 



Proof of Lemma 14.11 We expand the left hand side of equation (|4.4|) in terms of u{p, I) and 
rearrange the result. To this end we define the matrix-valued function X'-' : C — >■ c'""^^^'""^^ 
to be the (n — 1) x (71 — 1) submatrix of 

A A 
A A 

whose (1, 1) entry is the {I + l,r + j) entry. Then 

n 

G(p) = E'^('''')detX'^(p). (4.6) 



By Notation 12.71 and equation 1)4. 611 . the left hand side of equation ()4.4|l is equal to 



C(j^._i)/2(p)detX'-' - C(j^'_i)/2(p)^2jfcdetX' 

\^j-Jj odd j:J'^ odd k — l 



1 3 



1 3 



\3'--Jj even 



j:J'- even 



H{p). 



Splitting the sums over k into k : Jk is odd and k : Jk is even and rearranging inside the 
parentheses, we evaluate the square bracket to 



j:Jj odd y k-.j'i^ odd fc:J^ even / 

+ Y -cj^./2(p)e"''' - Ec(j/ _i)/2(p)lfcj +e"'''Ecj//2(p)Afc^ detX 

r-Jj even y k:j'^ odd k-.j'^ even / 

Making the change of variables fc 1— > r defined by 

Jfe is odd if and only if B r = {j'j, ~ l)/2, in which case k = , 
j'f; is even if and only if 9 r = Jfe/2, in which case k — J~ , 

it is clear that each of the parentheses in expression (|4.7|l evaluates to ^1 j . Hence 

n n 

jeJ+ jeJ- '=1 J=i 



n (4.7) 



1 3 



□ 



Proof of Lemma 14.21 The t-transforms of the boundary functions are entire, as are the mono- 
mials Cj, hence the sum of products of a f-transform and monomials Cj is also entire. By 

□ 



equation (|2.16|) this establishes that expressions 1)4. 5|) are entire functions of p. 



5 Alternative characterisations 

In this section we discuss sufficient conditions for well-posedness of initial-boundary value 
problems and present a proof of Theorem 11.51 These topics are unified by the arguments and 
notation used. 
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5.1 Sufficient conditions for well-posedness 



Throughout Section 15.11 we assume the boundary conditions are non-Robin. This simplifies 
the PDE characteristic matrix greatly, leading to corresponding simplifications in the argu- 
ments presented below. Nevertheless, we identify suprising counterexamples to the qualitative 
hypothesis 'highly coupled boundary conditions lead to well-posed problems whose solutions 
may be expressed using series.' 

We give the condition whose eflects are of interest. 

Condition 5.1. For A, a boundary coefficient matrix specifying non-Robin boundary condi- 
tions, we define 

C — \{j : oikj,Pkj 7^ for some k}\, the number of boundary conditions that couple the 
ends of the space interval, and 

^ = |{i ■ = for ^}|) number of right-handed boundary functions, whose 
corresponding column in ^4 is 0. 

Let a = ±i and let A be such that 

if n is even and a = 
ii±i if n is odd and a = i I + C. 
if n is odd and a — —i } 

We investigate the effect of Condition 15.11 upon the behaviour of the ratio 

J7m(p) 



ApDE (p) 



(5.1) 



in the limit p — >■ oo from within D. The PDE characteristic determinant is an exponential 
polynomial, a sum of terms of the form 



Z{p)e 



where Z is some monomial and Y C {0, l,...,n— 1}. As the problem may be ill-posed r]m is 
defined as in Lemma [3.11 a sum of terms of the form 



X(p)e-"'^«e^" / e-"""" qT{x)Ax (5.2) 
Jo 

where X is some monomial, qr G C°° [0, 1], F C {0, 1, . . . , n — 1} and z G {0, 1, . . . , n — 1} \ Y . 
Fix j G {1, 2, . . . , n} and let p £ Dj. Then the modulus of 

g->E„er"''p (5.3) 

is uniquely maximised for the index set 



Y 



{i - 1, i, • • • , J - 2 -I- f } n even, 

{ j - 1, i, . . . , j - 2 + i (71 + Im(a))} n odd. 



By Condition l5.1l ApnR (p) has a term given by that exponential multiplied by some monomial 
coefficient, Zj{p). That term dominates all other terms in Apde (p) but it also dominates all 
terms in r]m{p)- Hence the ratio (|5.ip is bounded in Dj for each j G {1, 2, . . . , n} and decaying 
as p — >■ oo from within Dj . 

If it were possible to guarantee that Zj 7^ then it would be proven that Condition 15.11 
is sufficient for well-posedness. Unfortunately this is not the the following example 

shows. 

Example 5.2. Let 

/I -1 0\ 
A = 1-10 
\0 1 2/ 

and consider the problem 11(3, ^4, i, 0, go). Then 



Di C |pGC:0<argp< II 
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and 

/ -C2{p){e-'P-1) -ci(p)(e-'^^l) -co(p)(e-^^ + 2) \ 

A{p) = -C2(p)(e--'' - 1) ~ujc,{p){e-^-<' - 1) -u^\o{p)(e-^^' + 2) . 

V-C2(p)(e-"^'''-l) -Uj\i{p){e-'^''P ~l) -iucoip)ie-'"'^ + 2)J 

We calculate 

ApDE (p) = (oj - oj')c2(p)ci(p)co(p) [9 + (2 - 2)(e''' + e'^^ + e'"'") 

+ (1 - 4){e-"' + e-'"" + c"'"'")] , 

in this case, as /3i 2 + /32 1 + /33 = 0, the coefficients of e"^^'' cancel for each j, 

= 3(0; - a;2)c2(p)ci(p)co(p) [s - (e-'" + e''^^ + fi-""'")] , 

2 

ri3{p) = {uj'^ - uj)c2ip)ci{p)co(p)'^uj^Tioj^ p){e" '' - e'''" - e""" '' + !). 

Fix J > 0. Consider a sequence, (pj)jgN, defined by pj = Rjc''^^^^ where 
00 

Rj ^ U ((1 ' ^) - 5, (1 - ^) ^Trm + 5) (5.4) 

m=0 

is a strictly increasing sequence of positive real numbers with limit 00 chosen such that pj G Di 
and. The ratio ri3{pj) / ApoE (pj) evaluates to 

^qrip,) - ujqTiujp,)e-''^'-''^f^ + cj^ric^'' Pj)e'^"^+''^P^ + 0(1) 



3^g-i(i-c^)Pj + 1) 4- o{e~'^o(^-'^')/^^) 

The denominator is 0(1) but, by condition (15. 4p . is bounded away from 0. The terms in 
the numerator all approach infinity at different rates, depending upon qr- Hence the ratio is 
unbounded and, by Theorem ll.il the problem is ill-posed. 

Indeed, third order initial-boundary value problems with pseudo-periodic boundary condi- 
tions are ill-posed if and only if 

a = i and /3i 2 + /32 1 + /Js = or 

,1 1 1 ^ 
a ^ -I and — + — + — ^ 0. 

pi 2 P2 1 P3 

A combinatorial necessary and sufficient condition for Zj 7^ in odd order problems is pre- 
sented as Condition 3.22 of [TJ] but is omitted here due to its technicality however we do 
improve upon that condition; see Remark 16.21 No further third order examples are known 
which obey Condition 15. II but are ill-posed. 

Condition 3.22 of 14 may be adapted to even problems by setting k = n/2 — R. The 
pseudo-periodic problems of second and fourth order are ill-posed if and only if 

n = 2 and = /3ii +/32 0, 

n = 4 and = /3i 3/32 2 + /32 2/33 1 + /^s i/34o + /34o/3i 3 + 2(/3i a/^s 1 + ^22l3io)- 
For example, the problem 11(4, A, ±i, h, go) with boundary coefficient matrix 



A : 



A 


1 

















M 








1 


-1 


























1 


1 








\o 

















1 


-1/ 



is ill-posed. 
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Remark 5.3. The essential difference between tlie odd and even cases presented above is 
that for odd order problems the well-posedness criteria depend upon the direction coefficient 
whereas for even order problems they do not. This means it is possible to construct examples of 
odd order problems that are well-posed but whose solutions cannot be represented by a series 
using Theorem 11.41 Indeed the problem Il{3, A,i, h, qo), with boundary coefficient matrix 
given by 

A = 

\o 1 2^ 

is well-posed but is ill-posed in the opposite direction. This is the issue mentioned in Re- 
mark 3.3 of [B]. 

Remark 5.4. There are classes of examples for which Zj 7^ is guaranteed. Indeed, Condi- 
tion [ST] is precisely the necessary and sufficient condition for well-posedness of problems with 
simple boundary conditions proved in [12j . 

Remark 5.5. There exist problems 11 for which 11' is ill-posed but for which 

^^^^^ ^ as p — ^ 00 from within , 




ApDE (p) 



for all j £ J"*" or from within E~ for all j £ J~ . This is a property of the (^j, dependent 
upon which column of A is replaced with the transformed data, not of the sectors in which the 
decay or blow-up occurs. In this case it is possible to deform contours over the corresponding 
hence one of the terms 

^ipx-ip"t (p) ^ f ^ip(x-i)-ip^t Cj (p) . 

aE+ ^+ ApDE (p) JsE- ^ ApDE (p) 



evaluates to zero in equation H1.8|l but the other does not. 

Remark 5.6. It is a conjecture that Condition 15.11 together with Condition 3.22 of [14] (as 
modified above to include n even) are necessary as well as sufficient for well-posedness of 
problems with non-Robin boundary conditions. Any counterexample must satisfy several 
strong symmetry conditions that appear to be mutually exclusive. Indeed for a problem, 
which fails Condition 15.11 or which satisfies Condition 15.11 but for which Zj = 0, to be well- 
posed several monomial coefficients X from equation (|5.2I) must be identically zero. 
Remark 5.7. We give a condition equivalent to Condition 15.11 for Robin type boundary condi- 
tions. Indeed, we define 

Bi — \{j £ J~ : 3 k,r for which P^j^O and ^ is a pivot}], 

B2 ~ IJ I and 

-B3 — \{j € : 3 fc,r for which a^j 7^ and /3fcr is a pivot}]. 
Then the condition is 

if n is even and a = 
^ if n is odd and a = i > ^ B2 + -Bs- 
^^^-ji if n is odd and a = —i ) 

5.2 Series representations for n even 

Proof of Theorem 11.51 By Theorem 11.11 the well-posedness of Il{n, A,i, h, qo) and the argu- 
ments of Section ISTTI for each j £ {1,2, . . . ,n} there exists some Fmax C {0, 1, . . . , n — 1} such 
that 

1. the term 

appears in Apde with ^y„,,,x 7^ a polynomial, 
ApDE (p) = O ('i2'y_(p)|e''"(''S«e^max - 
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2. for all y C {0, 1, . . . , n - 1}, z e {0, 1, . . . , n - 1} \ y for which 



is a term in rjjc, for some k, with Xy z 7^ a polynomial such that 
Xy . (p)e-"' ^ ^^^^^ QT (c^V) = o (1 (p) le'""^'' 



as /9 — > 00 from within Dj . 
Hence, for all such Y , z, 



Im 



i<f) 



< 



(5.5) 



for all X e (0, 1) and for all (p £ (0, 7r/n). 

If Il{n, A, —i,h,qo) is ill-posed then there exist Y, z satisfying the conditions above, x £ 
(0, 1) and (/) G (ji /n,2n /n) such that 



Af = Im 



> 



Define Fmax = {w - y ; y G Fmax} and 



Y\j{z} Im(e"*a;") > 0, 
y Im(e*'*a;^) < 0. 



Then, as n is even. 



Im 



hence there exists some x G (0, 1) such that 



> iV > 0, 



Im 



>0, 



which contradicts inequality (|5.5p . The argument is identical in the other direction, switching 
the intervals in which S lies. □ 



6 PDE discrete spectrum 

In this section we investigate the PDE discrete spectrum, the set of zeros of an exponential 
polynomial. We use the definitions, results and arguments presented in [11) . 

Lemma 6.1. The PDE characteristic determinant and PDE discrete spectrum have the fol- 
lowing properties: 

1. APDE (p) = (-1)"-'APDE M). 

2. /et y C {0, 1, . . . , n — 1}, y = 1 mod n : y £ Y} . Let Zy and Zyi he the polynomial 



coefficients of exp (—ip^ 



and exp {~ip'}2y^y/ respectively, in Apde . Then 
Zy[p) = {-ir-^Zy,(up). 

3. either Apde is a polynomial or the PDE discrete spectrum is asymptotically distributed 
in finite-width semi-strips each parallel to the outward normal to a side of a polygon 
with order of rotational symmetry a multiple of n. Further, the radial distribution of 
the zeros within each strip is asymptotically inversely proportional to the length of the 
corresponding side. 
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Proof, (i) The identity 

Akjii^p) = Ak+ij{p) 

follows directly from the definition H1.6|l of the PDE characteristic matrix. A composition with 
the cyclic permutation of order n in the definition of the determinant yields the result. 

(ii) By definition there exist a collection of index sets y C ^{0, 1, . . . , n — f } and polynomial 
coefficients Zy{p) such that 

Apde(p)= ^ ^y(p)e-'^^«e^'"". 
Yey 

By part (i), 

Define the collection y' ^ {{y + I mod n : y e Y} : Y e y}. Then 

J2 ZY{p)e-"'^y<^^''' = (-1)"-^ J2 ZY'ip)e-"''^y<^^"^\ 
Yey Y'ey 

Equating coefficients of exp {~ip'}2yeY '^^) yields y — y' and the result follows. 

(iii) The result follows from part (ii) and Theorem 8 of [TT]. □ 



An immediate corollary of Lemma 16. f I is that the PDE discrete spectrum has no finite 
accumulation point and is separated by some e > 0. 

Remark 6.2. A corollary of (ii) is that Zy = if and only if Zy' = 0. This means it is 
only necessary to check Zj ^ for a particular j in conjunction with Condition [STT] to ensure 
well-posedness. This permits a simplification of the general Condition 3.22 of |14) . 
ft is possible to strengthen part (iii) of Lemma |6. II in certain cases. 

Theorem 6.3. Let n ^ 3 be odd and let A be such that Apde is riot a polynomial. If n ^ 7 we 
additionally require that Condition \5. 1\ holds and the relevant coefficients, Zj, are all nonzero. 
Then the PDE discrete spectrum must lie asymptotically on rays instead of semi-strips. 

Proof. Assume n ^ 7 and the additional conditions hold. If 

'-^]--y 

then, by part (ii) of Lemma f6. II {0, I, . . . , (n — 3) /2} G 3^ hence the indicator diagram of Apde 
has subset the convex hull of 



S = <^ a;-- ^ cjf : r G {0, 1, . . . , n - 1} 

I yeY 

If F = {f , 2, . . . , (n + I)/2} G y then the indicator diagram contains the regular 2n-gon that 
forms the convex hull ofSUl — sisGS*}. We show that the indicator diagram is precisely 
the convex hull of S or of {—s : s G S} and that there are no points X^^gy '^"i i^o^ ^' ^ 
on the boundary of the indicator diagram other than at the vertices. 

Excepting rotations of Y and Y , which all correspond to vertices, the sets Y' £ y whose 
corresponding exponent has greatest modulus 

veY' 



are rotations and reflections of 



Fi = {l,2,...,(n-3)/2} or 

y2 = {f,2,...,(n-3)/2,(n + I)/2}. 
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However, the minimum modulus of the boundary of the indicator diagram is greater than or 
equal to 



(n-3)/2 



(n-l)/2 



E E 

y=0 y=l 



(n-l)/2N 



{n-3)/2 

+ E ^ 

y = l 





(n-3)/2 




(n-3)/2 




> 


E 


= Si > 


(n+l)/2 , a 


= S2 




y^l 




y = l 





hence any point corresponding to Y' is interior to the indicator diagram. It is easy to check 
that this also holds if n = 3 or n = 5. 

As there can only be two colinear exponents lying on any side of the indicator diagram, the 
argument in Sections 1-7 of [TT] may be simplified considerably and yield the stronger condition 
that the zeros of the exponential polynomial lie asymptotically on a ray, a semi-strip of zero 
width. The arguments of Sections 8-9 applied to this result complete the proof. □ 

Remark 6.4. Theorem 16.31 does not hold for n even. Indeed, 




n/2 

E^ 



(n-2)/2 

E ' 



hence if 

{0, 1, . . . , (n - 2)/2}, {1,2,..., n/2}, {1, 2, . . . , (n - 2)/2} G 3^ (6.1) 

then, by part 2. of Lemma [6. II there are three colinear exponents on each side of the indicator 
diagram. Condition (|6.1|l does not represent a pathological counterexample; it is satisfied by 
most pseudoperiodic, including all quasiperiodic, boundary conditions. 
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